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1. Introduction 

The sigma term and the strangeness content of the nucleon are phenomenologically relevant 
quantities which, however, cannot be directly measured in experiment. 

They are of interest, because they link the pion-nucleon and kaon-nucleon scattering lengths 
to the hadron mass spectrum, to the quark mass ratio m s /m ud , where m ud = {m u +m d ) /2 denotes the 
isospin averaged first generation quark mass, and - last but not least - to the strangeness content of 
the nucleon which, in turn, has a say on the importance of quantum fluctuations. 

That the sigma term cannot be directly 1 measured in experiment follows from the definition 

C nN ={N(p)\m u uu + m d dd\N(p)) = ^ (1.1) 

o k n = (N{p)\m ud (uu + dd)\N(p)) = m ud (1.2) 

am ud 

since, in nature, we cannot change the quark masses. Here we have specified the simplification 
that emerges in the isospin limit where the up and down quarks assume a common mass m ud . 
By contrast, on the lattice we can vary the quark masses, and this opens a unique opportunity to 
determine the pion-nucleon sigma term from lattice QCD datasets with several pion masses. 

Provided the lattice simulations include a strange quark (which ours do), one can also study 
the kaon-nucleon sigma term and the ss-nucleon sigma term 

1 1 . r 1 dM N dM N i 

Okn = -{m ud +m s ){N\uu + ss\N) = -{m ud + m s )\ -- h — \ (1.3) 

2 2 12 dm ud am s ) 

Oun = 2m s (N\ss\N) = 2m s ^ (1.4) 

am s 



where in practice it is common to trade (1.4) for the strangeness content of the nucleon 



2{N(p)\(ss)(0)\N(p)) 
yN (N(p)\(uu + dd)(0)\N(p)) 1 ' ' 

to which it relates vmy^m s /m ud = Os S n/Gkn- Also the sigma terms are linearly dependent in the 
isospin symmetric case (which we will assume in the following, unless stated otherwise), since 
GnN I m ud + GssN I m s = ^&kn / {p-ud +m s ). The sigma terms have the dimension of a mass (without 
showing any scheme or scale dependence), while the strangeness content is a pure number. 

On the lattice there are two main strategies to determine the sigma terms. One option is to stick 



with the definitions in ( |l.l[jl.4[ ), which then leads one to evaluate qq in a nucleon in and out state. 
This is technically involved (and noisy), due to quark line disconnected contributions [jj]]. The 
second option is to use the Hellmann-Feynman theorem (in a form adapted to field theory [Q]) and 
to determine the sigma terms from the variation of the nucleon mass as a function of the up/down 
or strange quark mass. In the following, we chose the second option, albeit in the version where 
one trades the dependence on m ud and m s for one in M\ and M? = 2M\—M\, as this choice avoids 



Measurable quantities like the nN and KN scattering lengths can be linked to (1.1, 1.2) and (1.3) by means of 



formulas from XPT, but this should not be mistaken as a direct measurement of the sigma term. 
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Figure 1 : Overview of our simulation points in terms of M n and Mg s = (2M% — M 2 . ) 1 1 2 . The lattice spacings 
for the three values of j3 are a~0.124fm (j3 = 3.3), a~0.083fm Q3 =3.57), and a~0.065fm Q3 =3.7), 
respectively, and the physical point is marked with a cross. Error bars are statistical only. 



renormalization issues. In other words, we measure for various combinations of M\ and M? 
(and various lattice spacings a and box sizes L), interpolate the results with ansaetze which we will 
discuss below, and evaluate the slope of the interpolation function in the relevant direction at the 
physical mass point. In practice, this interpolation is actually an extrapolation in M\ (whereas the 
data have essentially the right M?.), and it is clear that outside the range where we have data the 
uncertainty on the derivative grows much faster than the uncertainty on the function itself. 

In the following we discuss the details of our ensembles in Sec. 2, and some key features of a 
novel functional ansatz for octet baryon masses derived from covariant baryon chiral perturbation 
theory (CBXPT) in Sec. 3. First experiences with this formula and, as a complement, with more 
traditional polynomial and rational ansaetze are reported in Sec. 4 and Sec. 5, respectively. The 
main goal is to give a reliable estimate which precision can be achieved with our current data, and 
this together with some preliminary numbers and some outlook is arranged in Sec. 6. 



2. Overview of our "6 stout" dataset and scale setting issues 

The dataset to be used is the one that has been generated for our study of the hadron spectrum 
in QCD [Oh, and the overall spirit of the analysis is the same one as in the fx/ fit paper [ffl. 



3 



Sigma term and strangeness content of the nucleon 



S. Diirr and A. Ramos 



The ensembles have been generated with two degenerate light quarks and a separate strange 
flavor, a scheme commonly referred to as A^y =2+1 QCD. They are based on a Symanzik action with 
tree-level coefficients for the gluons and a clover action with 6 levels of "stout" or "EXP" smearing 
applied to all links which enter the covariant derivative or the clover term. This combination was 
found to entail very good scaling properties in standard hadron observables []|]. Note, finally, that 
our multiply smeared fermion action is as local as the unsmeared clover quark in the usual sense 
of locality, i.e. D(x,y) = for ||jt— y\ \ > 1. Still, in the sense of gauge field locality it has a more 
extended range: 8D(x,x)/8U^(y) is non-zero for \\x— y\\ larger than 1 unit, albeit with a steep 
fall-off pattern - see the supplementary online material of [J3j] for details. 

Our ensembles include three lattice spacings (a ~ 0.065, 0.083, 0.1 24 fm) and the light quark 
mass m ud is varied such that the pion mass covers the range from 190 MeV to 670 MeV. By contrast 
the strange quark mass m s is almost the physical one; the combination M? = 2M\ — M\ is always 
in the vicinity of its physical value. The physical box size L is such that M n L is in the range of 
4 or larger; this limits finite-volume effects on our hadron masses to an amount smaller than the 
statistical fluctuations [|J]. An overview of these ensembles (to which we refer to as the "6 stout" 
or "6 EXP" data) in terms of M n and Mg s = (2Mf — M^) l l 2 is given in Fig. [I]. The smallness of our 
minimal pion mass M™ 11 ~ 190MeV bears the promise that the extrapolation to the physical point 
is a relatively mild one, thus entailing a controllable systematic error in the final result. 

All dimensionful quantities in the previous paragraph implicitly build on knowledge of the lat- 
tice spacing a. This brings us to the issue of how the scale is set which, as we shall see, receives an 
extra twist when derivatives w.r.t. the quark mass are taken, relative to the case when only spectral 
quantities at the physical mass point are calculated. The physical mass point is the point where any 
2 of the 3 ratios that one may form from aM n , aM^, aMx (with X being the particle through which 
the scale is set, we will consider X = Af,S,£2 below) take on their physical values. Apart from the 
quantity X to be used to set the scale, there is also a choice regarding the scale setting scheme. In [^] 
we used two such schemes, the "mass independent scale setting scheme" and the "ratio method". 
In the former case for any given j3 the measured aM^, aMx, aMx are interpolated by a smooth 
function of the bare quark masses am ud , am s and at the point where (aM % ) / (aMx ) , (aMx ) / [aMx ) 
assume their physical values the interpolated aMx is identified with a times the physical value of 
Mx ; this yields the lattice spacing a for all ensembles with a common coupling parameter j3 . In the 
latter case dimensionless ratios are formed on a per ensemble basis, and these ratios are interpolated 
with a smooth function of am ud , am s , and read off at the point where (aM % ) / (aMx ) , (oMk ) / (aMx ) 
assume their physical values. Effectively this means that the scale is set for each ensemble indi- 
vidually; the lattice spacing a depends on the combination (j8 ,am ud ,am s ). While these two scale 
setting schemes yield identical results (at the physical mass point) for all spectral quantities [g], 
there is a slight subtlety if sigma terms are evaluated via the Hellmann-Feynman theorem. 

For definiteness, let us consider the pion nucleon sigma term at a fixed physical value of m s , 
as a function of m ud . The way it is calculated on the lattice amounts to the factorization 



dM " „ ,„2 dM N 



OnN(m ud ) = m ud ~ M 



dm ud n dMl L a JpdgLm x dMl 



Ml dM N 



(2.1) 

latt/phys— pt 



where the apparent dependence on the scale setting channel X actually boils down to cut-off effects 
at the physical mass point (which then disappear in the continuum limit), while it seems there is an 
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ambiguity in the last bracket for which the continuum limit does not provide a remedy. The point 
is that the derivative in the last bracket makes reference of the true scale away from the physical 
mass point, and this is, as we have just seen, ambiguous. Perhaps it is easiest to discuss an extreme 
case first. If we set the scale through the nucleon, then the pion nucleon sigma term with the "ratio 
method" is - by definition - zero, i.e. a^ mt = 0. And if we set the scale through the omega, the 
pion nucleon sigma term with the "ratio method" evaluates to the same amount as the difference 
between the pion nucleon sigma term and the pion omega sigma term in the "mass independent 
scale setting scheme" would, i.e. a^, rat = a^ sss - a^ sss . 

Taking a look at (LI, Q| ) reveals the origin of this apparent discrepancy. The matrix elements 
refer exclusively to the physical mass point, and can be evaluated without any ambiguity (at least 
with a Ginsparg- Wilson type action). The not-so-innocent part is the derivative w.r.t. the quark mass 
taken in the last equalities. Here it is (implicitly) assumed that the lattice spacing does not change 
as the quark mass is varied, and this is why we must assume - for consistency reasons - the "mass 
independent scale setting scheme" when evaluating sigma terms, on the lattice, via the Hellmann- 
Feynman theorem. In other words, the requirement to attribute a common lattice spacing a to all 
ensembles with one coupling /3 comes about through the specific transcription of the observable of 
interest; in general there is no restriction on the scale setting scheme. Still, it holds true that the 
QCD j3 -function depends, for asymptotic coupling (i.e. for small enough lattice spacing), only on 
the number of active flavors, not on their masses. Accordingly, the "mass independent scale setting 
scheme" is always applicable, since it stipulates a property at finite lattice spacing which, due to 
asymptotic freedom in QCD, must hold at arbitrarily small lattice spacing. 



3. Octet masses in CBXPT 

Given the discussion in the previous sections, it is clear that a controlled determination of the 
pion nucleon sigma term relies on accurate measurements of the nucleon mass at various pion 
and kaon mass points, and a smooth ansatz = M^(M^ t ,M} s ) which is valid in the entire regime 
from the heaviest datapoint included in the fit down to the physical pion mass (and eventually down 
to the chiral limit, if one wishes to determine chiral low-energy constants). 

In the present contribution we focus on o k n at the physical mass point. For this purpose it is 
fully sufficient to come up with an analytic (i.e. polynomial or rational) expression for M^(M^,M?) 
and an analogous expression for Mx(M%,Mj s ) for the state X that is used to set the scale, since - in 
the interval between M^ hys and M™ ax of the dataset - QCD is an analytic function of the quark 
masses. Preliminary results from such an approach will be reported in Sec. 5 below. 

In a further perspective it is clear that we will not be able to resist the temptation of testing the 
host of predictions by chiral perturbation theory (XPT) of how different quantities relate to each 
other. In XPT the quantities connect through their behavior in the (2-flavor or 3-flavor) chiral limit, 
and this is the reason why, with this goal in mind, the ansatz must be good all the way down to zero 
quark mass. The A^-flavor chiral limit of QCD is dominated by the logarithmic singularity induced 
by the spontaneous breaking of the flavor SU (N/)a symmetry and the dynamics of the pertinent 
pseudo-Goldstone bosons. For observables built from pseudoscalar mesons the consequences have 
been cast into a valid form in the seminal papers by Gasser and Leutwyler [||, ^. 
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Table 1: Summary of the baryon-meson-baryon couplings gMB in terms of the low-energy constants D,F. 



For observables that include baryons the task is significantly more involved, and we attempt to 
give a short explanation why there is a plethora of different versions of chiral perturbation theory for 
baryons. Ultimately, they all rest on Weinberg's power counting theorem but the basic difficulty 
is that the quadratic divergence in self-energies, brought about by a meson loop attached to a given 
line, cannot be traded for a logarithmic divergence, as is the case in the dimensionally regulated 
setup of meson XPT. The original approach taken by Gasser, Sainio and Svarc |g] ("OBXPT") 
is to live with this fact and to pay a certain price in terms of a slowly converging series away 



from the chiral limit. The heavy baryon approach by Jenkins and Manohar [ 1C] ("HBXPT") treats 



the nucleon as a non-relativistic particle, and in this approach it was shown that the A can be 



introduced as an explicit degree of freedom and that this improves convergence [11, 12]. The 



covariant approach taken by Becher and Leutwyler [13] ("CBXPT") aims at establishing relativistic 
covariance again, via infrared regularization. To the best of our knowledge, in this framework other 
octet members have not been included as explicit degrees of freedom yet. 

In short one can say that the HBXPT approach is very successful for those observables where 
the physics is completely dominated by light (i.e. physical or lighter) pions which are treated as 
the pseudo-Goldstone mesons in QCD with 2 light flavors. Extensive work in the nineties showed 
that applications of HBXPT to baryon observables that depend on the strange quark flavor lead, 
in general, to large (i.e. unnatural) cancellations between leading order and next-to-leading order 



contributions (see e.g. [14, 15, 16]). During the first decade of the new millennium, this insight 
was painfully re-discovered by the lattice community, as several collaborations attempted chiral 
extrapolations (m uc j — > at fixed m s ) of 2+1 flavor baryon data, using only leading-order HBXPT 
formulas. This triggered a search for "better" (i.e. more practical) recipes, for instance 

• finite-range regulator extensions of HBXPT to include or model higher order effects (see e.g. 
Amherst group [|l5|] and Adelaide group j jl7| , |l8|]) 

• taming higher order effects by adding a host of decuplet contributions in HBXPT/NRSSE 



(see e.g. LHP Collaboration []19p and others [gfj]) 

resort to polynomial/rational fit functions in M 2 oc m u d (e.g. us [^, ||] and other groups) 

1 /2 

resort to polynomial/rational fit functions in M K oc mj d (see [[l^, pi]]). 



modify integration contour in CBXPT approach (cf. Dorati-Gail-Hemmert [22]) 



but a fair review of these developments is clearly beyond the scope of this contribution. 

In the following, we are going to explore the suitability of one such formula, which belongs to 
the last point in this list. It has been worked out by one of us (TH) and involves the function 



H(X 2 ) = -^r { J\-^-~ arccos (^-) + i^-log (^) \ (3.1) 
v ; 4n 2 [V AM 2 \2M J 4M & \M 2 ) J 



6 



Sigma term and strangeness content of the nucleon 



S. Diirr and A. Ramos 



where X will be a meson mass. Furthermore, it builds on SU (3) relations among the various 
baryon-meson-baryon couplings gMB, as is evident from Tab. |[ The set (D,F) relates to the more 
common pair (§,£a) through £ = F/D and gA = D + F. The two main attractive features, from 
our point of view, are that this formula describes how the complete set of baryon octet states varies 
as a function of the meson masses, and that there is a chance that it remains adequate up to higher 
pion masses than is usual in the HBXPT approach (say up to M K ~ 400 — 500 Me V rather than 
M % ~ 200-300 MeV). In full glory it reads 

M N =M - 2(2b + b D + b F )Ml - 2{b Q + b D - b F )Ml 

+^H(Ml) + ^H(M 2 K ) + ^H[M\) +A8 n Ml + A8, s Mi (3.2) 
F £ Fr F n 

M A = Mo - 4(b + b D /3)M 2 K - 2(b + 4£ D /3)M? 

+^H(Ml) + 8 -^H(M 2 ) + ^H(M 2 ) + 4M4 + 4feM* (3-3) 

^71 **rj 

M z = M -4(b Q + b D )M 2 n - 2M4 

+ 8 -§H{Ml) + 8 -§H{M 2 K ) + ^H(M 2 ) +45X + 45^ (3.4) 

M z = M - 2{2b + b D - b F )M 2 K - 2(b + b D + b F )M% 

+^H(M 2 K ) + 8 -^H(M 2 K ) + 8 -^H(M 2 ) + A8 K M A K + 4^4 (3-5) 

with all meson-baryon couplings gMB parametrized by only two constants, as indicated in Tab. 
In total this gives a parametrization with 8 unknowns: Mo, bo, bp, b F , gA, 8 n , 8s S - Note that the 
coefficients in front of the M\ and M% contributions are common to all octet members. To the order 
we are working at only two of the three pseudo-Goldstone boson masses are linearly independent; 
they are connected through the Gell-Mann-Okubo relation 3M^ = 4M|- — M 2 . Accordingly, it 
makes sense to consider (M 2 ,M 2 S ) the basic mass coordinates. 



4. First experiences with CBXPT extrapolations 



Before applying formula (3.2-3.5) to our data, let us make a few practical comments. First, it 
is important to notice that this formula builds on SU (3) chiral symmetry. Accordingly, Mo denotes 
the (common) mass of the baryon octet in the 3-flavor chiral limit. Phenomenology suggests a value 



in the range M ~ 770 MeV, although with a large error margin [23]. Clearly, this is an opportunity 
for the lattice to come up with a considerably more precise determination. Next, the coefficients 
D and F are reasonably well known from phenomenology. It makes sense to fix the combination 



gA = D + F = 1.2694(28) [ |24] ] to its value at the physical nucleon mass, as the difference to the 
value in the chiral limit is expected to be small. For the ratio ^ = F /D the situation is less clear; 
there are two preferred scenarios in the literature, ^ = 2/3 and ^ ~ 0.5. 

A few comments are in order regarding the pseudoscalar decay constants F k ,Fk,F^. First, let 
us clarify that we use the "Bernese" normalization where F^ hys = 92.2 MeV. The more relevant 
point is that, whatever choice is made for it is not supposed to destroy the SU (3) chiral 

symmetry in the 3-flavor chiral limit. Accordingly, pinning F k ,Fk,F 1] down at their phenomeno- 
logical values is not an option. We see three legitimate choices: (z) use a joint 3-flavor chiral 
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CBXPT: IVU=0.402GeV,x =[8.39,0.25], a^O.OSSGeV 
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Figure 2: Result of one particular "snapshot fit" of the ansatz (3.2-3.5) to our "6 stout" dataset with a mass 
cut M n < 410MeV. For display purposes the data have been shifted, by means of (4.1), to [M?] p h ys (top) 
and [M^jphys (bottom), and only the remaining dependence on M\ (top) or M? (bottom) is shown. 
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limit value Fq, (ii) use an SU (3) chiral formula (e.g. the one at NLO with Fo, 14,1,5) (Hi) use the 
measured data (perhaps with an SU (3) compatible interpolation, in any case this requires Za). 

In the following we present the result of one particular fit based on choice (/). This fit addresses 
the full baryon octet and yields, in principle, eight sigma terms ((T^jy, o&jy and ditto for A,£,E). 
Due to its preliminary nature, and because the comparison with the analytic approach is not ready 
for other channels, we will just quote Ojin, OssN- Still, one should keep in mind that these numbers 
are preliminary, and the assessment of systematic uncertainties is not yet finalized. 

Having made these cautionary remarks, we are in a position to present in Fig. the result of a 



fit of the ansatz (3.2-3.5) to our data. The parameters Mq, bo, bo, bp, 8 n , 8 Ss , Fq have been adjusted 
by the fit (to reasonable values), while (£,£a) = (2/3, 1.2694) was held fixed. All together, this 
is a fit with 7 free parameters to 40 datapoints with an uncorrected % 2 = 8.39. The resulting 
uncorrected # 2 /d.o.f. ~ 0.25 seems plausible, since the four octet masses in each ensemble being 
highly correlated will lead to an underestimate of the true # 2 /d.o.f. by about a factor four. 

Note that this is the result of one "snapshot tit", i.e. with a specific choice of the fitting window 
[fmin,fmax] for each state, with one pion mass cut (here M % <410MeV), and so on. This particular 
fit yields g k n = 55(10) stat MeV, where the quoted error is only statistical, and y N ~ 0.16. To get a 
trustworthy estimate of the systematic uncertainty, one should consider reasonable variations over 
the fitting range, the pion mass cut, the scaling behavior of cut-off terms, and the functional ansatz 
for the dependence on (M 2 ,M} S ), as will be briefly discussed in Sec. 6 below. 

A technical point worth mentioning is the shift recipe applied to show the result of the fit. As 
is clear from the discussion, the fit spans, for each baryon octet member, a two-dimensional surface 
above the (M 2 ,M 2 S ) coordinates depicted in Fig. |l[ However, because it is difficult to graphically 
display how 40 datapoints behave relative to 4 surfaces, we resort to two one-dimensional plots. 
We "shift" the data, along the surface established in the fit, to the physical value of M 2 S and plot the 
remaining dependence on M 2 . In other words what is plotted in the first panel of Fig. ^ is [Q] 

data(M 2 , 2M 2 -M\) - fit(M 2 K , 2M 2 -M 2 K ) + fit(M 2 , [2M 2 -M 2 K ] phys ) (4. 1) 

and a similar shift is applied, in the second panel, to bring all datapoints to a common value of M 2 
and depict the remaining dependence on M 2 S . It goes without saying that the recipe ( fi~T| ) does not 
affect the fit itself, it just helps to display the result. 

5. First experiences with analytic extrapolations 

A complete analysis of a phenomenological variable at the physical mass point must include 
a variation over the fitting ansatz that is used to interpolate or extrapolate the data. To this aim we 
include polynomial and rational ansaetze into our analysis, in the same way as we did in [Q]. 

This analytic approach builds on the fact that an expansion of QCD Green's functions about 
the physical point (m^ ys ,m? hys ) is completely regular. Therefore it makes sense to define dimen- 
sionless expansion parameters A^ ~ (m ul j — m^ s )/A, A gs ~ (m s — m^ hys )/A, and to express the 
measured octet mass as a function of these parameters. At this point details matter. We prefer 

I aMji ^ {M* hys \ 2 ( aU- ss \ 2 (Mf y \ 2 
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POLY: o =0.053(1 4)GeV, RATIO: G jiN =0.044(6)GeV 
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Figure 3: Result of one particular "snapshot fit" with the polynomial and rational ansatz (5.3,5.4), where 
the scale is set through the N. The same shift procedure (4. 1) has been applied to the data as in Fig. 2. 



with X being the scale setting state over the alternative choice 



A' = ( 



aM n 



KaMx 



M 
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A' = ( 



\aMx 



\2 /M? hys \ 
) ~\~M^) 



(5.2) 



since in ( p.l| ) all dependence on the simulated am ut i, am s is in the very first numerator, while in 



( p.2| ) the numerator and the denominator of the first terms depend on the quark masses. 

With such small mass parameters in hand, and bearing in mind that the range of simulated M\ 
is much larger than the range of simulated M| (cf. Fig. fil), we consider the polynomial ansatz 



a-M) 



phys 



l + c 1 A 7t + c 2 A 2 7[ + c 3 A s - s 



(5.3) 



where c\ ,02,03 and the isolated a (per coupling j8, cf. the discussion in Sec. 2) are the fit parameters. 
Obviously, with this choice one sets the scale through the nucleon mass, which means that the 
quantity that is effectively calculated is o n ^/M N at the physical mass point. Likewise 



a-M] 



phys 



1 - diA K - d 2 A\ 



■ d^A- s 



(5-4) 



is a rational ansatz with similar characteristics. With any such fit in hand, one proceeds along the 
lines of (2.1 ) and evaluates the change of the fit function under a change of (aM K ) 2 , (aMg s ) I . With 
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POLY: M<=0.402GeV,x =[8.16,0.34], a reN =0.052GeV 
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Figure 4: Result of one particular "snapshot fit" with the global polynomial ansatz, where the scale is set 
through the fi, withM^ < 410MeV. The same shift procedure (4.1) has been applied to the data as in Fig. 2. 
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a simple ansatz like (5.3) or (5.4) this can even be done analytically, giving 



dM h 



m ud 



dm ud 
dM N 



M 



im s 
m s Gs sN 



mud GkN 



phys 



phys 



phys 



■2M[ 



dM N 
9Mk 



Mi 



M N 



phys 



mi 



s dMl 



phys 



C3 



2c 3 /ci 



(5.5) 

(5.6) 
(5.7) 



In the same spirit, one may consider a global polynomial or rational ansatz which fits the entire 
baryon octet at once. In this case setting the scale through an octet member would introduce some 
asymmetry, and it seems advisable to include one more fit to use X = SI for this purpose. 

With the same cautionary remarks applicable as in Sec. 4, we show some "snapshot fits" to 
the ansaetze ( p3| ) and ( gg ) in Fig.|. They yield c kN = 53(14) stat MeV and o nN = 44(6) stat MeV, 
respectively, whith the quoted errors being statistical only. The result of a joint "snapshot fit" to 
the full octet is presented in Fig. |j. Here, the scale is set through the SI, and the same pion mass cut 
M K <410MeV is used as in Fig.^. Both g k n = 52(10) stat MeV and ~ 0.13 are well compatible 
with what was found in the CBXPT approach. In fact, with either fit the strangeness content is 
well consistent with zero. Going back to Fig. [T] and the bottom panel of Fig. || or Fig. |] it appears 
that this is linked to our strange masses being slightly above the physical target value. Likely, with 
simulated values of 2M\—M\ covering the entire range between (say) 0.3 GeV 2 and 0.65 GeV 2 
[such that the physical value of 0.47 GeV 2 would be the center of this range] one would have a 
better chance of determining the slope. Note that for the other octet members, the situaton looks 
better, since their slope in the ss direction seems to be larger. 



6. Preliminary results and outlook 

Given the discussion in the previous two sections, it is clear how we should proceed to come up 
with a valid determination of the nucleon sigma terms (or equivalently of o n N, yjy) and analogous 
sigma terms for the baryon octet states. 

We have implemented several functional ansaetze to describe the dependence of the baryon 
state on (M 2 ,M?). Both the CBXPT ansatz tested in Sec. 4 and the family of polynomial and 
rational ansaetze tested in Sec. 5 yield reasonable results for the sigma terms (and in the former 
case also for the low-energy parameters). With such an interpolation in hand, one may compute 
the derivatives with respect to M 2 , M? , at the physical mass point. In some cases this is a simple 
function of the fitted parameters. In other cases it proves more convenient to evaluate the derivatives 
with respect to the measured M % , Mk and to convert via 



O nN (m u d) ^Ml 



dM N 



dM\ 



Ms S fixed 



Mi 



dM N {M K ,M K ) 2 dM N {M K ,M K ) dM\ 



dMl 



+ M: 



dM 2 K 



dMl 



Mg S fixed 



M K dM N {M n ,M K ) , Ml dM N (M Kl M K ) 



dM K 



+ 



4M K dM K 



(6.1) 



dM, 



N 



dMl 



Mx fixed 



2M, 



2 dM N (M n ,M K ) dM 2 K 



dMl 



dMl 



Mx fixed 



Mj dM N (M Kl M K ) 
2M K dM K 



(6.2) 
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where again everything is evaluated at the physical mass point. 

What remains to be done is the systematic variation of all ansaetze over the original fitting 
window, over the pion mass cut, and over the scaling pattern of potential discretization effects. Of 
course, all of this should be done on 0(1000) bootstrap samples to assess the statistical uncertainty, 
but this part is standard. This kind of machinery was used in Ref. [Q] to give a reliable assessment 
of both the statistical and the systematic uncertainty of the observable of interest {/k/ fit)- 

For the pion-nucleon sigma term in MeV units our fits usually yield values in the lower fifties 
with typically about ten MeV statistical error. With hindsight we thus anticipate that the final 
result will be in the range of g k n ~ 50(10) (10) MeV, where in each slot only one digit is meant 
to be significant. Regarding a^N or the situation is less convincing. With the "6 stout" dataset 
depicted in Fig. [j] we obtain, with each ansatz, large statistical errors and non-negligible spreads 
among the ansaetze. Currently, a value like — 0.1(2) (1) seems appropriate, which would not 
even tell whether there is a non-zero strangeness content at all. Our current understanding suggests 
that, in order to obtain a substantially more precise value, one would have to add simulation points 
with significantly smaller strange quark mass than we have right now. 
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